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Abstract 



We study the conditions under which an overdamped regime can be at- 
tained in the dynamic evolution of a quantum field configuration. Using a 
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real-time formulation of finite temperature field theory, we compute the effec- 
tive evolution equation of a scalar field configuration, quadratically interacting 
with a given set of other scalar fields. We then show that, in the overdamped 
regime, the dissipative kernel in the field equation of motion is closely re- 
lated to the shear viscosity coefficient, as computed in scalar field theory at 
finite temperature. The effective dynamics is equivalent to a time-dependent 
Ginzburg-Landau description of the approach to equilibrium in phenomeno- 
logical theories of phase transitions. Applications of our results, including a 
recently proposed inflationary scenario called "warm inflation" , are discussed. 

PACS number(s): 98.80 Cq, 05.70.Ln, ll.10.Wx 

I. INTRODUCTION 

Kinetic equations describe the time evolution of a certain chosen set of physical variables. 
The choice of physical variables in principle is arbitrary, but often in practice is governed 
by the measurement of interest. Typical examples are the order parameter of a complex 
system or the coordinate of a Brownian particle in a heat reservoir. The kinetic approach is 
usually implemented through a proper separation of the microscopic equations of motion of 
the chosen physical variables into regular and random parts. An averaging over the random 
part then generates the effective partition function for the regular part. This averaging is 
often referred to as a coarse-graining. 

One typical application of the kinetic approach is when the physical variables of interest 
possess energy in relative excess or deficiency to the rest of a large system. Kinetic theory 
then describes the approach to equilibrium of the chosen physical variables, as for example 
in the kinetics of phase transitions or in Brownian motion. In the former case, the system is 
able to release energy to the environment due to some change in its internal state. Provided 
the environment is disproportionately large relative to the system, the process is irreversible. 
For a continuous transition, the focus of the present work, this process of equilibration can be 
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described by the monotonic change of an appropriate order parameter, which is the chosen 
physical variable. Many systems are known to relax in this manner. Phenomenologically, 
they are successfully described by the time- dependent Ginzburg-Landau theory [[[]. 

Here we are interested in examining under which circumstances physical variables whose 
microscopic dynamics is second order in time, as for example the Higgs order parameter 
of spontaneous symmetry breaking, may have a dynamics which is effectively first order in 
time as in Ginzburg-Landau phenomenological models. 

Qualitatively it is not difficult to argue the plausibility of this standard picture for the 
Higgs symmetry breaking scenario. A single variable, the Higgs order parameter, is modeled 
to control the release of energy to all the modes that couple to it. By basic notions of 
equipartition, one anticipates that some portion of the order parameter's energy will flow 
irreversibly to any given mode. Provided the Higgs order parameter couples to a sufficient 
number of modes, the motion of the order parameter will be overdamped. 

In particle physics models, Higgs symmetry breaking is accompanied by mass generation. 
Thus the natural couplings for the Higgs field <fi to bosonic fields Xi is 4> 2 xl > gauge fields 
is 4 >2 A lfl Ai f j i and fermionic fields ipi is 0^^. For a microscopic realization of time dependent 
Ginzburg-Landau theory for the Higgs scalar order parameter in a particle physics setting, 
these are the most obvious types of couplings to investigate. In this paper we will examine 
the case of purely bosonic couplings in the "symmetry restored" regime. That is, we will 
study the relaxation of an order parameter which is initially away from the only minimum 
of the free energy density describing the system. Much of the formalism required for this 
has already been done in but we will extend that calculation to the overdamped regime. 
In an upcoming paper, we plan to study the symmetry broken case. 

To our knowledge, this paper is the first study of overdamping in quantum field theory 
with realistic couplings between system and environment, as inspired by particle physics. 
Overdamping has been studied in quantum mechanical reaction rate theory for a particle 
escaping from a metastable state (for a review please see |3]]). This is sometimes referred 
to as the Kramer's problem, with the overdamped limit also called the Smoluchowski limit. 
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Quantum mechanical models describing this problem are commonly of the system-heat bath 
type. Microscopic quantum mechanical models have been constructed along these lines, in 
which the particle (system) is coupled to a set of otherwise free harmonic oscillators (heat 
bath). Such microscopic system-heat bath models are often referred to as Caldeira-Leggett 
(CL) models. In many cases they have been exactly solved ||. The overdamped limit has 
been derived in these models for the case where the coupling is linear with respect to the 
oscillator variables but arbitrary with respect to the particle variable 

A Caldeira-Leggett type model has also been formulated for the case where the system is 
a self interacting scalar quantum field coupled linearly to a set of otherwise free fields and the 
overdamped limit has been obtained ||. This model does provide a microscopic quantum 
mechanical realization of time-dependent Ginzburg-Landau dynamics in scalar quantum field 
theory. However, since the couplings between system and environment variables are linear, 
it should be considered as a first step toward more realistic treatments. More importantly, 
the calculational method used in cannot be extended to the case when the system variable 
couples quadratically to other fields. 

Although the analysis of overdamping in this paper has general applicability, it was 
motivated by the warm inflation scenario of the early universe [§,0. In |7]] it was realized 
that the standard Higgs symmetry breaking scenario, when put into a cosmological setting, 
provides suitable conditions for the universe to enter a de Sitter expansion phase and then 
smoothly exit into a radiation dominated phase. The overdamped motion of the order 
parameter in this scenario may sustain the vacuum energy sufficiently long for de Sitter 
expansion to solve the horizon and flatness problems. Simultaneously, the relaxational 
kinetics of the order parameter can maintain the temperature of the universe and permit 
a smooth exit from the de Sitter phase into the radiation dominated phase. Finally, the 
thermal fluctuations of the order parameter provide the initial seeds of density perturbations, 
which in addition could be scale free under specified conditions . An elementary analysis 
of this scenario, based on Friedmann cosmology for general realizations of order parameter 
kinematics, indicated that if the universe's temperature does not fall too much during de 
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Sitter expansion, then the cosmological expansion factor from the de Sitter phase should be 
of order the lower bound set by observation ||. Although this is not a tight constraint of 
this scenario, it is a natural one. An analysis of COBE data motivated by this expectation 
did indicate a slight preference for a small super-Hubble suppression scale, which could be 
interpreted as arising from a de Sitter expansion with duration near its lower bound [JTD 



Furthermore, the overdamped limit required by warm inflation, when expressed in different 
terms, was noted || to be an adiabatic limit, for which known methods from dissipative 



quantum field theory P, [illJl4j1 are presumed valid. These facts provide further motivation 
to seek a microscopic model of the scenario, which is the goal of the present work. 

The calculational methods used here, based on Schwinger's close-time path formalism, 
were developed in [[J. There are several other works in the literature that apply this formal- 
ism to a variety of different situations. [See, for example, the works of Refs. [|12|-|T6|1 .] The 
new feature of the present paper is to shift focus to a kinematic regime dominated by strong 
dissipation, in order to establish under which conditions this regime leads to overdamped 
motion. This approach will allow us to have a unique understanding of the microphysi- 
cal origin of such dynamical behavior, which is in general invoked phenomenologically in 
applications ranging from condensed matter physics to inflationary cosmology. 

The paper is organized as follows. In Sec. II our model of interacting bosons is presented 
and the effective action is computed perturbatively for a homogeneous time dependent back- 
ground field configuration <p{t). In Sec. Ill the effective Langevin-like equation of motion 
is obtained for <fi in the symmetry-restored phase. In Sec. IV the overdamped limit of this 
equation of motion is derived and regions of validity are given. In Sec. V the results of 
the previous sections, which are for Minkowski space, are extrapolated into a cosmological 
setting and a preliminary examination is made of the warm inflation scenario. In Sec. VI 
concluding remarks are given. Two Appendixes are included to clarify a few technical details, 
like the evaluation of the imaginary part of the self-energies and to stress the importance 
of taking fully-dressed field propagators to properly describe dissipation in the adiabatic 
approximation for the field configuration. 



II. MODEL OF INTERACTING BOSONIC FIELDS 



A. The Effective Action 



Let us consider the following model of a scalar field <fi in interaction with N scalar fields 



Xj- 



where 



1 N rl 1 

m = g (<V>) 2 - V[<f>] + £ {-(d, Xj ) 2 - V[ Xj }\ - VUt, Xj] 



v[<t>] = + j/ 



(2.1) 



(2.2) 



and 



v\xA = fx? + %x) 



v-^xh = ji 9 h\) 



(2.3) 



(2.4) 



For the most part, we will consider all coupling constants positive: A, fj and #J > 0. Writing 



— ► (p + 1] in (|2.1|) , where (p is a background field configuration and 77 are small fluctuations 
around (p, we obtain the expression for the 1-loop effective action T[<p], valid to second order 
in the fluctuations, by performing the functional (Gaussian) integrations in 77 and \f 



N 



I>] = S[<p] + liTr In [□ + V"{ip)\ + \i £ Tr In □ + + ^ 



2 vr/J 2 

where S[</?] = ftfxCfaO], V"(<p) = 



(2.5) 



and 



A? 



^Tr In [□ + V"{cp)} + ^ Tr In □ + /i? + g*<p 
2 j=i 



In y Dr? n D Xj exp {-^ [□ + V"(<p)] 77 - [ n + + 3> 



(2.6) 



Neglecting contributions to Q2.5| ) which are independent of we can expand the loga- 
rithms in (|2.5| ) in powers of <p, obtaining, in the graphic representation: 
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Xj Xj 

n^i-ioop = X3 + XX + EiLi I X) + XX 1 + ° (A3) + ' (2 - 7) 



where we have identified the propagators in the internal lines. External lines are Xip 2 /2 for 
the 0-graphs and g 2 <p 2 for the x-graphs. 

B. Single-Particle Excitations and Dissipation: Dressing the Propagators 

Before presenting our derivation of the effective nonequilibrium equation of motion for ip, 
we contrast our approach with earlier works in the literature. We closely follow the method 
of Ref. || in the derivation of the evolution equation for (p. In particular, it was shown in |fj] 
that for slowly changing fields, dissipative terms vanish if they are computed perturbatively 
with bare propagators. There are several issues related to this result. Boyanovsky et al. 



in Ref. |13| argue, in the context of a toy model, that dissipative effects cannot be studied 
within perturbation theory: perturbation theory breaks down before dissipative effects can 
be observed. This shows that dissipation is a nonperturbative effect in quantum field theory. 
In it was shown that dissipative terms can be derived once a consistent "dressing" of 
propagators is used. This is an explicit way of considering the effect of quasi-particles (or 
single-particle states) in the evolution of the system, described by ip, in interaction with 
a thermal bath which represents fluctuations of <fi and of others fields to which it may be 
coupled. 

It seems reasonable to expect that dissipation effects are closely related to the effect of 
collisions which dress the field propagators. Take for example the case of a bare propagator 
expressed in terms of the spectral density Po{p), where there is a one to one correspondence 
between the energy and the momentum of a given state. This completely neglects the 
spreading of possible energy states due to interactions. In a full "dressing" of propagators, 
this is accounted for through the introduction of a lifetime (decay width) for single-particle 
states, such that the full (dressed) spectral density p{p) is smeared out. In particular, particle 
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lifetimes are crucial in the study of relaxation time-scales in quantum many-body theory 



HI- 

Also, the reason why we can get dissipation within our approach can be traced back to 
the very way that transport coefficients are derived in quantum field theory. As we will 
show later, the assumption of a slowly moving field is consistent with overdamping in a 
strong dissipative environment, justifying the adiabatic approximation we adopted. In this 
regime, there is a close relation between the dissipation we compute and the shear viscosity 
computed from the Kubo formula [|i9H22|. As explained in |20|j2l[| , diagrams contributing to 
the shear viscosity have near on-shell singularities for free bare propagators. Full resummed 
propagators regulate these singularities through an explicit thermal lifetime of single particle 
excitations. Analogous singularities are exhibited by our expressions for dissipation terms 
if bare propagators are used. Additional issues concerning the relation of our dissipation 
terms with the shear viscosity will be discussed in the following two sections. 



C. Self-Energies and Dressed Propagators 



From the above discussion, we rewrite the Lagrangian density in fl2.1|) as 



2 V ' -*j-r if* +2 



+ £ {^te) 2 - \ {fi + E*) x) - f <?x) - f ix) + \^x]\ , (2.8) 

where and S Xj . are the self-energies for the and Xj fields, respectively. This way we can 
work with full (dressed) propagators for the 4> and Xj fields (note the implicit resummation 
of diagrams involved in this operation) and at the same time keep consistency by considering 
|y0 4 — |£<^ 2 and also ^Xj ~ \^xjX 2 j as interaction terms. This method has already been 
adopted before in many different contexts [see, for example [@,^,^]]. In terms of the self- 
energies the field propagators are written as 

1 1 



q 2 — w? + ie q 2 — m 2 — £(g) + ie 



(2.9) 
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For both and Xj fields, a finite lifetime of single particle excitations, given in terms of 
the imaginary part of the self- energies, first appear at the two-loop order. We thus restrict, 
for simplicity, the evaluation of and S x , up to the two-loop level. Diagrammatically, the 
self-energies are given by 



X X 



X 



0( )0 x()X x()x 



X 

+ ( ^ ) + ( £ \ + higher loop terms ( 2 -10) 





and 



X 



X 9 9 X 



X 



+ ( x ) + f x ) + higher loop terms. (2.11) 




X 



The setting sun (non-local) diagrams in ( |2.10| ) and (|2.11|) (the two last terms in ( 2.10 ) 



and ( |2.11| )) contribute imaginary terms to the self-energies, from which we can write the 
decay widths T^, T Xj , for the and Xj fields, respectively, in terms of the on-shell expressions 
(£/ = ImE): 

r, (?) = (2 . 12) 

and 

r Xj (g)= 4^ ' (2 - 13) 

where u;^ ( x ) is given by the solution of u; 2 = q 2 + m 2 + Re£(q, u). 

Explicit expressions for T(q) in the A0 4 model have been obtained in |2l[ and [22]. We 



follow |2T[ to compute and T Xj . A straightforward extension of the computation can be 
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applied to our model of interacting 0— \j fields. Some of the details are shown in Appendix A, 
where we evaluate the imaginary contribution coming from the mixed setting sun diagrams 
in and S Xj . [the last diagrams in ( 2.10 ) and ( p.llj )]. Even though in general there are 
no simple way of expressing the results, if we adopt the zero space momentum (|q| = 0) 
approximation for the imaginary part of the self-energies, we can find simple approximated 
expressions for both (|2.12| ) and ( p,13|) , respectively, given at finite temperature (/3 = l/T) 
by (for m T ~ 0(^(T))) 

X 2 T 2 



r*(?)l 



Li 2 1 



-(3m,T 



and 

r x ,(g)| 



- (1+2^* 



f 2 T 2 



N 



Li 2 (l - 



-0m T \ 



) -Li 2 



(2.14) 



E/(°>WCD) 2 8 7r 3 cu x .(q) 



f 1 + \^. X] 



Li 2 1 - e 



9^ 



2 5 tt 3 ^ X3 (q 



Li 2 1 



Li-; 



'1 _ e -/3M 3 (T)- 



1 



-/3m.T 



[2.15) 



In the above expressions, itlt and /ij(T) are the thermal masses for and Xji respectively 
S^ )Xj = 1 for rriT = Hj(T) and 8^ Xj = otherwise. Li 2 (z) is the dilogarithm function^. 

This approximation for the decay widths, in terms of the zero space-momentum expres- 
sion for the imaginary part of the self-energies, is common to computations of transport 
coefficients and contrast densities in field theory [11,1^,20]. However, Wang, Heinz and 
Zhang |22] showed that this approximation may lead to errors in the calculation of the con- 
trast density in the A0 4 model. In fact, the expressions for ImE can be fast changing for 
some momentum range and values of the masses. For example, in Fig. 1 we plot the value of 
(the on-shell) ImS(g), obtained numerically, as a function of the momentum, normalized by 



1 We follow the convention in Ref. |24|] for the definition of the dilogarithm function: Li 2 (z) = 
— f* jzjdt. Some useful approximations for Li 2 (z) are Li 2 (z) ~ ^- + [hi(-z) — M z + C(-2 2 )) f° r 
z < 1, and Li 2 (z) ^ - \ ln 2 (z) + 0{l/z), for z > 1. 
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its |q| = expression (for fj <C gf). Even though ImE(g) can depart considerable from its 
|q| = value, we will show later that, for a range of small thermal masses, this approxima- 
tion results in a small error (< 10%) in the expression for the dissipation coefficient, when 
compared with the computation using the complete |q| ^ expressions for ImS(g). 

In the analysis presented in the next sections, it will also be sufficient to use the lead- 
ing order high temperature expressions for the finite temperature effective (renormalized) 
masses, uit and /ij(T), appearing in (|2.12|) - (|2.15|) , (obtained from the 1-loop diagrams in 
( p.lO| ) and (|2.11|) , respectively), given byQ 



\rp2 N rp2 

m 2 T = m* + Re£,(m T ) T> Z m m 2 + ^ + £ 2% (2.16) 



and 



tf(T) = $ + BeZMT)) ^+ fj ^- + 9^ ■ (2.17) 



D. Real-Time Full Field Propagators 

In order to obtain the evolution equation for the field configuration tp, we use the real- 



time Schwinger's closed-time path (CTP) formalism |25|. In the CTP formalism the time 
integration is along a contour c from — oo to +oo and then back to — oo. For reviews please 
see, for example, refs. f26|-|2"8|. 



In the CTP formalism the field propagators are given by (with analogous expressions 
for G v 



'XjJ- 



2r 



The divergences in ( |27i0| ) and ( ^TTT]) , as in the effective action, can be dealt with by the usual 



introduction of the appropriate renormalization counterterms in the initial Lagrangian, for the 
masses, coupling constants and the wave-function. In particular, we note that the imaginary terms 
in the self-energies expressions, coming from the setting-sun diagrams, are finite, m, fi, g, f and A 
in ( 2.16| ) and (2.17) and in our later results are to be interpreted as the corrected and not as bare 



quantities. 
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+ x ) 

G^ (x, x ) 
(x, x ) 

^6 ( X 1 X ) 



i(T + <f)(x)<f)(x')) 
i(T-<f>(x)<f>(x')) 

mx')^)) 
mx)^')) , 



(2.18) 



where T + and T_ indicate chronological and anti-chronological ordering, respectively. G~^ + 
is the usual physical (causal) propagator. The other three propagators come as a conse- 
quence of the time contour and are considered as auxiliary (unphysical) propagators. The 
expressions for G r ^ l {x : x') in terms of its momentum-space Fourier transforms are given by 



G,p(x,x') 



(2vr) 3 



G+ + (q,t-if) G+-(d,t-t 
G£+(q,t-f) G-f^t-t* 



) 



(2.19) 



where 



G+ + (q,t-t') 
G-f{^t-t') 
G+-(q,t-t') 
G7 + (q,t-t') 



G>(q,t- t')6(t - + G<(q,t - - t) 

G>(q £ - 0«(f " *) + G<(q, t - 0«(t - f ) 

Gj(q,t-0 
^(qt-f) 



(2.20) 



In terms of the decay width T^, the expression for the full dressed propagators at finite 
temperature where obtained in 0], from which we have 



G>( q ,t - 0= 7T - {[1 + ™H - ^)] e- 4 ^-^^*-*') + n(^ + zT )e^+ ir ^-*')} 



2a;0 

G<(q,t -*')=£> (q,t'-t) 



(2.21) 



-i 



where ra(c<j) = (e^ w — l) is the Bose distribution and uj = uj(q) is the particle's energy, or 



dispersion relation, kty(q) = yq 2 + For G x ., w x (q) = Wq 2 + /i](T) 



III. DISSIPATION IN THE ADIABATIC REGIME 
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A. The Effective Equation of Motion 



With fields in the forward and backward segments of the CTP time contour identified 
as an d 4>-,X-i respectively, the classical action can be written as 

S[<j>,x] = I d A x{£{<P + , X+ ] -£[0-,X-]} , (3.1) 

The evaluation of the effective action at real time can be done exactly as in M. There 
are also a number of other works using Schwinger's closed-time path formalism to obtain 
the real-time effective action for field configurations. [See. e.g., Refs. [0-0.] Here we will 
concentrate on the evaluation of the effective equation of motion in the strong dissipative 
regime. In the evaluation of the effective action there appear several imaginary terms, 
once the x fields and the fluctuations around the (p background are integrated out. These 
imaginary terms can be interpreted as coming from functional integrations over Gaussian 
stochastic fields, as can be visualized by introducing the new field variables: 

<-Pc = ~ (0+ + 0-) , <^A = 0+ - 0- • (3-2) 



In terms of these new variables the equation of motion is obtained by P,|14 



5<f A 



\ifiA=0 



where £j are stochastic fields, related to each distinct dissipative kernel appearing in (|3.3j ). 



At 1-loop order, the leading contributions to the dissipative terms in the equation of 
motion come from the diagrams: 

>CX +E MXX1 (3-4) 

Xj 

The explicit expression corresponding to these terms appearing in the effective equation 
of motion, Eq. (|3.3| ), is (as obtained in for a similar case) 
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N 



(p c (x) \ d 4 x'<p 2 c (x') 



-Im 



G 



++ 



G 



A 2 f d 3 q n (l + n ) " g) f d 3 q n Xj (l + n Xj ) 



+C A 2 



8^7 (2vr)3^ (q)r ^ (q) ^ 2 -7 ( 27 r)3^(q)r Xj .(q) 



4 X Xj 



Xj. 



d 3 q J A 



2V 



(2vr) 



— Im 

2 



G++(q,t-t')J +2^4lm[Gj+(q,i-0 



(3.5) 



where in the lhs of the above equation, we used the compact notation 



C++ 



-l 2 

x,x' 



d 3 k 
(2^)3 



exp [ik.(x — x')] 



d 3 q 



G++(q,t-t')G++(q-k ! t-t'), (3.6) 



(2tt) 3 *>* 

with G ++ (q, t — t') obtained from ( |2.2U| ) and ( |2.21| ). In the rhs of (|3.5|), we have taken 
the limit of homogeneous fields, for details see the Appendix B. We have also made use 
of the approximation for slowly moving fields: <f 2 (t') — <$.(t) ~ 2ip c (t)ip c {t)(t' — t). In the 
next section we show that this approximation is consistent with strong dissipation. After 
performing the time integration and retaining the leading terms in the coupling constants, 
we obtain the result given in ( |3.5| ). The last term, proportional to (p 3 , will correspond to 
the finite temperature correction to the quartic (f) self-interaction (see Appendix B) . 

The final equation of motion, at leading order in the coupling constants, at high tem- 
peratures (/ij(T), ttit <C T) and in the adiabatic limit, can then be written as 



<p c + m 2 T ip c {t) + ^¥>c(*) + W 2 c{t)<Pc(t) = ¥> c (t)fi(0 > 



(3.7) 



where ttit is given by fl2.16j ), At is the temperature-dependent effective (renormalized) quar- 
tic coupling constant^: 



3 The terms linear in the temperature come from the two- vertex diagrams in (|3.4|). The apparent 



instability from these terms for high T is only an artifact of the loop expansion. As shown in [29] 
for the (j) model, once higher order corrections are accounted for, At is always positive even in the 
T — > oo limit. Using full dressed propagators we are automatically taking into account these higher 
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T 



A T ~ A - — { — : 

2 l87rm T 

N t T 

+ 0(A 3 ,^/,AV,^ 6 ), 



+ o 



hi 



4ttT 



+ 7 



+ o 



T 



+ 



(3.8) 



In ( |3.7| ), £i is a stochastic field associated with the imaginary terms in the effective 
action coming from the real-time evaluation of the diagrams ( |3.4| ). Its two-point correlation 
function is given by [0 



<6(*)a(x')> 



A : 



N 



Re 



G\ 



3=1 



C++ 



2 

x,x' 



(3.9) 



Note that since we are considering homogeneous field configurations, £i is a space uncor- 
rected stochastic field, but it is colored (time dependent) and Gaussian distributed, with 
probability distribution given by (iV\ is a normalization constant) 



= JV^exp { ~ I d 4 xd 4 x'^{x) 



A 5 



-Re 



G\ 



G 



Xj 



&G0 



(3.10) 



As shown in 0, the dissipative coefficient in ( p.7|) , written explicitly in Eq. (|3.12 ) below, 
and the noise correlation function Eq. (|3.9| ) (in the homogeneous limit), are related by a 
fluctuation-dissipation expression valid within our approximations (1-loop order at A 2 ,g| 
and for T/uj < 1 , Y/T < 1): 

^Jd 4 x'UxMx')). (3.11) 



In it was also shown that as T — > oo, r<^ )Xj — > oo, and the integrand in ( |3.9| ) becomes 
sharply peaked at |t — t'| ~ 0. In this limit, we can obtain an approximate Markovian limit 
for 



order corrections, through the appearance of thermal masses in fl3.8| ). However, in the multi- field 
case there is the possibility of vacuum instability due to the <fi couplings to the Xj fields. This 
appears as a constraint in our estimates below. 
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We can read the dissipation coefficient 771, which appears in Eq. ( |3.7|) , from (|3.5|)P1, 



„ = —8 [ d3q + n <t>) , V 9 J-B [ n ^ 1 + n ^ , n ( A 2 — n 4 ^-) 

Vl 8 P J (2vr)3 ^(q)r (q) + ^2 P J (2vr)3 ^ (q)r% (q) + [ ^ 9j ^ J ■ 

(3.12) 

For the model we are interested in, with Lagrangian density given by ( |2.1| ), with a large 
number of x fields coupled to 0, and for /j -C g 2 and A < gj, we can use the obtained 
expressions for and T X] , to show that 3> T Xj . . Since the dissipation coefficient, Eq. 
Q3.12j ), goes as 1/r, T x . will give the dominant contribution to rjx. An explicit expression 
for 771, can be obtained by using the |q| = approximation for ImE^g) and ImE Xi (g), or, 
equivalently, Eqs. ( |2.14| ) and ( |2.15| ) for and T Xj , respectively, in Eq. ( |3.12| ). At the high 



temperature limit, T ^> mjn,fij(T) and for ~ 0(fij(T)), we then obtain the following 
approximate expression for r)\ (using Li2(z) ~ tt 2 /6, for z Cl) 



r»m^(T) 96 J A 2 /2T\ 



KrriT/ 



+ S/l + B^l!-!,^)] 1 "^)}' <3 ' 13) 
In order to test the validity of the above approximate expression for 7^, we have computed 
it numerically. The two expressions are shown in Fig. 2, for fj <C g 2 , A < g,,-, and iV = 25, 
where, for simplicity, we have also considered [ij = /i and gj = g for all Xj fields (mj ~ 
5/ij(T)). We see that the above approximation for rjx fits reasonably well the full expression 
for the dissipation coefficient in the high temperature region, having a < 10% discrepancy 
for m T /T < 0.4. 



4 In H an extra contribution to the cj) decay rate coming from the (f> — X interaction was left out. 
Here we give the correct expressions for I?,*, T x and for the dissipation. 
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B. Dissipation Coefficient and Shear Viscosity 



It is interesting to note the close relation of the above expression for the dissipation 
coefficient with that obtained for the shear viscosity evaluated, e.g., from a Kubo formula 
MM: 



Vsheax 



ft 

rf 3 x / dt / 

oo J —oo 



dt / ([7r w (0),7r fcI (x,t , )]> 



(3.14) 



where itki = (5^5/ — g5*5/)T^-, with T*- the space components of the energy-momentum 
tensor. In our case, with Lagrangian given by ( j2.1[ ), 



Tf, 



dC dC 



(3.15) 



In order to compute the shear viscosity in (|3.14 ) to lowest order, we must evaluate the 
diagrams ( p.4|) , which, as shown in |20|,21|, have near on-shell singularities coming from the 
product of (bare) propagators. These singularities are softened once explicit lifetimes for 
excitations are included through dressed propagators. Taking this into account, we obtain 
the following expression for the shear viscosity r] shea _ r (in analogy with the evaluation of //shear 
in the A0 4 single field case), 



//shear 



T>m TlAtj (T) (3 f d?k 



30 J (2tt) 3 ' 



7^(l + n^) + n Xj (l + n x .] 



uSi r v . 

X] X] 



(3.16) 



Compare the above expression with ( |3.12 ). The evaluation of (p. 16 ) leads to the standard 
result for the shear viscosity being proportional to T 3 and inversely proportional to the 



coupling constants. However, Eq. (|3.16|) , as shown by Jeon in |21| , does not represent 
the unique contribution to r/ s hear at this order of coupling constants. Due to the near on- 
shell singularities and the way they are regulated by the thermal width, there is an entire 
class of diagrams, called ladder diagrams (diagrams with insertions of loops between the 
two propagators in (|3.4j) ), contributing to r/shear at the same order. By using a formal 
resummation of vertices, Jeon was able to perform the summation of the whole set of ladder 
diagrams in the simple A0 4 theory, showing that the true value of the shear viscosity is 
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about four times larger than the one loop result in the high temperature limit. Since our 
expression for the dissipation coefficient exhibits the same properties of r] shea _ r , we expect that 
these higher loop ladder diagrams will also give a significant contribution to the value of rji in 
( p,12| ). However, as we are dealing with the more complicated situation of several interacting 
fields, we will not attempt here to evaluate these contributions. From the example of the 
shear viscosity calculation in the single field case, these ladder contributions will only add 
to the one- loop result for the dissipation coefficient, not changing qualitatively our results. 
Thus, Eq. ( p,12| ) represents, at least, a lower bound for the dissipation, applicable in the 
strong dissipation regime, as we will show next. 



IV. ADIABATIC APPROXIMATION AND STRONG DISSIPATION 

We now investigate the validity and limits of applicability of our main approximations, 
in particular the adiabatic approximation. In order to arrive at the expression for the 
dissipation, Eq. ( |3.12|) , and to write the equation of motion for ip c as in Eq. Q3.7|), we 



assumed that the field (f c changes adiabatically [see (p.5| )]: 

tp 2 c (t') — tp 2 c (t) ~ 2{t' — t)ip c (t)ip c (t) + higher time derivative terms . (4.1) 
This approximation for the field configuration has recently been the focus of some attention 



in the literature | lq| . The authors in [15], working with soft field modes set by a coarse 
graining scale k c , showed that the adiabatic approximation breaks down once the field con- 
figurations (soft modes) oscillate with the same time scale as the dissipative kernels (with 
time scale given by ~ k" 1 ). However, here we work in a very different context. We are 
mainly concerned with the overdamped motion of the homogeneous field configuration tp c , 
i.e., when its oscillatory motion is suppressed. Therefore, the dynamic time-scale for (p c 
must be much larger than the typical collision time-scale (~ T _1 ). Note that this is a much 
stronger condition than the simple requirement that the field should change slowly in time, 



with time scale set by the frequency uj(k) = y k 2 + m\. Thus, we must examine when the 
condition 



— >r^ 1 (4.2) 

is satisfied. 

We choose T as the smallest of the two thermal decay widths r^, T x ., as it will set the 
largest time-scale for collisions for the system in interaction with the thermal bath. 

Note that in the evaluation of the dissipation coefficient in ( |3.5| ), the leading contribution 
to the first time derivative of ip c is of order As discussed earlier in connection with 

the shear viscosity coefficient, the dependence of the dissipation coefficient on the decay 
width r comes from using it as the regulator of on-shell singularities present in (|3.4|) at first 
order in the time derivative. In Appendix B we present an argument justifying the need 
of regulating with the decay width and also compute the next order contribution in the 
adiabatic approximation, showing the consistency of the results. 

Since the stronger the dissipation the more efficient the adiabatic approximation, the 
parameter range where ( [4.2] ) is valid leads naturally to the regime where tp c undergoes 
overdamped motion (in the sense of Eq. (|4.4|) below). If we consider the ensemble average 
of the equation of motion ( |3.7| ): 



5<p A 



1 1/9 A =0 



, (4.3) 



where (...) means average over the stochastic fields, then we define the overdamped regime 
when the (averaged) background configuration ip c satisfies 

7}i<pl<p c + m 2 T ip c + -jftp 3 c = . (4.4) 
o 

We also restrict our study to the high-temperature and ultra-relativistic region: T ~ 
|q| 3> rriT,Hj(T). We take the couplings gj, fj such that g 2 3> fj. Also, for simplicity, as 
before, we take all gj = g. At high temperatures we can then write for ( |2.16p and ( |2.17p 
(T 2 > 24m 2 /A, 12fi 2 /g 2 ) the expressions 

T 2 



(A + 2iV) - , (4.5) 

where is the number of fields coupled to <f), and 

19 



(4.6) 



A. Results for three different cases 

We will examine the condition for strong dissipation with overdamped motion for three 
particular choices of parameters, showing that there is a region of parameter space consistent 
with this regime. Using (|4.4j ), we can write the equivalent expression for (|4.2j) : 



m 2 T + f <^ 



Vi<Pc 



« r x . (4.7) 



In the estimates below, we evaluated both 771 [from Eq. (|3.12|) 1, and Y x (computed at 
|q| = T) numerically. The three cases analyzed are: 
Case 1: A ~ g 2 : In this case we obtain that 

r4~(l + 2JVy^, (4.8) 

Note that the last condition is written as a constraint for the positivity of Ay. With 
these values and for the case N = 25, we obtain the results shown in Fig. 3a, where we have 
plotted both sides of Eq. (|4.7| ). The region of parameters satisfying Eq. ( |4.7|) is given by 
the intersection of the region below the solid lines (the function T x ) with the region above 
the dashed line (\<p c /ip c \ computed for different values of ip c ). 

Case 2: A ~ g: As above, this is shown in Fig. 3b. The region satisfying Eq. ( |4.7|) is given 
again by the intersecting region below the solid line and above the dashed lines. 

In both Figs. 3a and 3b, the results are shown up to the value of uit satisfying the 
constraint for the positivity of Ay. 

Case 3: At ~ g 4 : This case follows a slightly different philosophy, of fixing the corrected 
coupling as opposed to the bare coupling. We have, 
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A« 9 * + 2^ = X( g ,N), (4.10) 

m 2 T ~(\(g,N) + 2Ng 2 )^, (4.11) 
with the additional constraint, 

\(g,N)<l. (4.12) 

The results for this case are shown in Fig. 3c, with the same interpretation as in cases 1 
and 2: the region satisfying Eq. (|4.7|) is given by the intersecting region below the solid line 
and above the dashed lines. The results are shown up to the value for ttit satisfying the 
condition (ETT2D . 



We note that the case A^ m g A is the one with the broadest range of validity in parameter 
space, as seen in Fig. 3c, followed by the case A « g, shown in Fig. 3b. For A = g 2 , the 
condition for adiabaticity is only possible for fairly large field amplitudes, which may be 
beyond the validity of a perturbative evaluation of the effective action. We will come back 
to this issue in the next section. In any case, we stress that there are several regimes where 
the adiabatic approximation is valid. 

In all cases, the smaller N the smaller the region of parameters that satisfies (JO). 



In particular, for iV < 2, we find no parameter range satisfying ( fOj) and therefore, the 
adiabatic approximation. This is consistent with the intuition that dissipation is caused by 
the decay of the <fi field into x fields and is more efficient the larger the number of decay 
channels available. We also obtain that ip c is always somewhat large ( > 2T) for the range 
of physical parameters satisfying (f4.2|) , for both cases analyzed, being even higher for case 
1. 

If in ( f4.7| ) we use instead of T Xj , the region of parameters improves considerably; since 
r^, ^> Y Xj for large N, it allows much smaller values of <f c /T. It should be recalled that 
determines the relaxation time scale for the <fi field. 

Finally, as discussed earlier, the expression we quoted for rji gives only a lower bound 
for the dissipation coefficient. As in the case studied by Jeon in [ffiH , higher loop ladder 
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diagrams can lead to a considerably higher value for rji. For several interacting fields, simple 
estimates show that these ladder diagrams scale at most as N. Therefore, they may well 
be of the same order as the leading order 1-loop contribution to the dissipation coefficient, 
given by the x _sec tor. We leave a more detailed analysis of the contributions coming from 
ladder diagrams to a future work. Additional contributions to the dissipation coefficient in 
( |3.12j ) only improve our estimates, enlarging the region of parameter space satisfying the 
adiabatic approximation; the ratio if c /T decreases, broadening the conditions under which 
the field undergoes overdamped motion (strong dissipative regime). 

It is worth mentioning that the <p ~ X coupling constant in Eq. (|2.1| ) can be negative and 
this also leads to interesting results. As an illustrative example, consider an even number of 
X-fields with the sign of the 4> — x coupling distributed so that 

2 k 

Vmt = £(-!)Wx? (4.13) 
and fj = f,j = l... 2k. In order that the potential be strictly positive, it requires 

\24 24 4 J K 1 

which for large N implies g 2 < f/6. In the alternating sign regime (ASR) the thermal 
masses are 

m 2 T » A T 2 (4.15) 

and 



/4^T 2 . (4.16) 
Following an analysis similar to above and for case 3 (At ~ g A ), we find a solution regime 



within the perturbative amplitude expansion, g<j) < X(f) < T, for g 2 < / 3 / 2 ln(2y 24//)/46 
and iV ~ l/g 4 - For example, these conditions are satisfied for / < 1.0, g 2 i$ 1/20. In this 
example A ~ g 2 , but this can be modified in several ways. In general, when the <fi — x 
couplings are distributed between positive and negative strengths, it controls the growth 
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of rriT due to the cancellation of thermal mass contributions from the x-fields. Restricting 
the magnitude of uit, in turn, increases the parameter regime and duration of overdamped 
motion. This example demonstrates another regime of overdamped motion in our model for 
small field amplitudes g 2 ^ 2 < Mt- 



B. Summing up the whole 1-loop series - The effective potential 



The fact that overdamping in (|3.7| ) for much of the parameter space demands large field 
amplitudes, at least within the approximation scheme used here, is a direct consequence 
of having a field dependent dissipation rj(<f) ~ ip 2 . Since in ( |2.7|) we are considering a 
perturbative expansion for the 1-loop effective action in the field amplitudes (that is, in 
powers of \<p 2 /2 and gfal), the need for large field amplitudes may place doubts on the 
validity of our calculations for a considerable portion of the parameter space. Below we 
address this issue in two different ways; first by comparing our results with an improved 
one-loop approximation and then by using the subcritical bubbles method JJIJ to test the 
validity of the effective potential for large-amplitude fluctuations. 

We start by computing the analog of ( f4.4p in the context of the whole one-loop approx- 
imation, i.e., when \<p 2 /2 and gfac i n are taken as part of field-dependent masses. 



For this, let us give an alternative computation of the evolution equation for (p c in terms of 
the tadpole method of Weinberg [|l3| , [3ll|3"4| : in the shifting of the scalar field, = <p c + 77, 



the requirement (if) = leads, at the one-loop order, to the evolution equation for cp c (for 
homogeneous fields) 

(p c + mVc + + ^c(v 2 ) + E#>c(x?> = , (4.17) 

j 

where (rj 2 ) and (xj) are given in terms of the coincidence limit of the (causal) two-point 
Green's functions G^ + (x, x') and G+ + (x, x'), respectively, which satisfy, in the fully dressed 
propagator form (see, e.g., Ringwald in 



G++(x, x') + / ^zE^x, z)G+ + (z, x') = i5(x, x') (4.18) 
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and 

□ + M • + G++(x, x') + J d'z^ X] (x, z)Gg(z, x') = i5(x, x') , (4.19) 
where, in ( f4.18|) and ( |4.19| ), E^(a;, x') and T, x .(x,x') are the (causal) self-energies for the 



<f) and Xj fields, respectively. By expressing r](x) and Xj{ x ) i n terms of mode functions, 
we can then evaluate the averages in ( [4.17| ). An explicit expression can be obtained in 
the approximation (equivalent to the adiabatic approximation) d^(y9 c )/c<^(<^ c ) <C 1 and 
c ^x(v 9 c)/^(v 9 c) "C 1, for which there is a WKB solution for the mode functions of the fields. 
In this paper, however, we will not carry out this calculation. A detailed study of this, in 
the context of an expanding background and along the proposals made in the next section, 
will be presented in a forthcoming paper. 

For now, we can present the result of this calculation, by using the simplest formula- 



tion proposed in [ll|, based on a relaxation-time approximation of the kinetic equation, 



for the calculation of the averages in ( |4.17| ). We can then show that the (ensemble aver- 



aged) evolution equation for <f c {t) can be expressed, in the quasi-adiabatic approximation 
(hydrodynamical regime of [[11]]), by 

^c + ^W+W^-O, (4.20) 

where V^ s ((p c ) = ^q^^ , is the field derivative of the 1-loop effective potential, 

where u;| = q 2 + + |<^ 2 and oj^. = q 2 + /x 2 (T) + g]^p 2 c are the field dependent frequencies, 
with masses given in terms of the thermal ones|, Eqs. ( |2.16|) and ( |2.17|) . Also, r/i in (|4.20|) 



5 Note that this will lead to the daisy corrected effective potential. In particular, once the thermal 
masses are being introduced in the derivative of the effective potential, expressed in terms of one- 
loop tadpole graphs in the Weinberg method, it is well known that this method leads to a consistent 



finite temperature effective potential |35|| , with daisy graphs incorporated. 
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is the same as in ( |3.12|) , but now with the masses replaced by the field dependent ones. 



In terms of Q4.20| ), in the overdamping approximation, the condition ( |4. 2| ) becomes 

< T . (4.22) 



Using Eq. ( f4.22| ), in the high temperature approximation for the fields, 
m ( f > (T)/T ) m X:j (T)/T -C 1, we can show that the results obtained earlier, in terms of the 
amplitude expansion for the effective action, for instance, the results expressed in Fig. 3 
(with rriT replaced with the field dependent mass m^(T)), remain approximately the same, 
for the cases where (p c ^ 2T. Thus, at least for these values of the field amplitude, higher 
order corrections do not add to the effective potential. In other words, at leading order in 
the high-temperature expansion, the field derivative of V e s can be just expressed as in Q4.4D, 

Kir ~ m T^c + At/6 <p 3 c - 

We can also address the issue of high- amplitude fluctuations by adopting a method 



suggested in Ref . |31[ , where it was applied to test the validity of the 1-loop approximation to 
the electroweak effective potential. We note that the results from this approach are entirely 
consistent with nonperturbative computations based on lattice gauge theories performed by 
K. Kajantie et al. [|32]] . 

The interactions of the field with a thermal environment will promote fluctuations 
around the perturbative vacuum. The subcritical bubbles method models these fluctuations 
as unstable spherically-symmetric configurations with a distribution of sizes and amplitudes. 
For details see Refs. [pip!]]. Using a distribution function for these configurations, it is 



possible to compute the RMS amplitude of the fluctuations [31 



(p{T) = VV>T-<¥>> T , (4.23) 



where (. . .) is the thermal average defined in Ref. [ST]]. For the effective potential of Eq. 



([4.4]), we obtain, 



<p 2 (T) ~ \m T T . (4.24) 
o 
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Since the perturbative approach for the computation of the effective potential relies on a 
saddle-point approximation to the partition function, it will only be valid for small-amplitude 
fluctuations about the perturbative vacuum. For potentials which exhibit spontaneous sym- 
metry breaking, it is customary to choose the maximum amplitude to be at the inflection 
point, f max ~ Pint- Here, since we have a potential with positive-definite curvature, we will 
conservatively assume that the perturbative expansion is valid for fluctuations dominated 
by the quadratic term of the effective potential, that is, for 

2 

^ iax ^ 12^ . (4.25) 

At 

The condition for the validity of the 1-loop approximation for the effective potential is 
then written as 

f(T) < v^Lx • (4.26) 

It is straightforward to apply this condition to the 3 cases analysed above. Since case 3 
is the one with a larger range of parameters satisfying the adiabatic condition, we use it as 
an illustration. From Eqs. (|4.10| ) and ( |4.11|) , we can write Eq. (|4.26|) , after dividing by T 2 , 



as the inequality 



lfig - N)] " 2> T^' (427) 



where, f(g, N) = 24m^/T 2 . This condition is easily satisfied for a large range of parameters. 
In particular, for A = 0.5, g = 0.3, N = 25, which are values inside the region of parameters 
allowed for overdamping shown in Fig. 3c for ip c > 2T, we obtain, (p ~ 0.3T and <^ max — 
16T, well within the range of validity of the small-amplitude approximation. We thus 
conclude that it is possible to attain the adiabatic limit of strong dissipation within the 
1-loop approximation scheme adopted here. 

V. APPLYING STRONG DISSIPATION TO WARM INFLATION 



The calculation in Sees. II-IV presented a microscopic quantum field theory model of 
strong dissipation in Minkowski spacetime. This section addresses the application of this 
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calculation to the cosmological warm inflation scenario. Although we will not present a 
detailed extension of our previous results to an expanding spacetime, we will argue that 
most of the modifications are quite straight-forward up to the requirements for the warm 
inflation scenario. 



A. Formulation 

Consider the standard Friedmann cosmology with Robertson- Walker metric 

dr 2 



ds 1 = dt l - R\t) 



+ r 2 d6 2 + r z sin 2 



(5.1) 



1 — kr 2 

We restrict our analysis to flat space, k = 0, and quasistatic de Sitter expansion, H = R/R pa 
const. For notational convenience, the origin of cosmic time is defined as the beginning of 
our treatment. For this metric, the minimally coupled Lagrangian for the model in Eq. (|2.1|) 
is 

" v d 3 xe 3Ht ((a o 0(x,t)) 2 - (e~ ffi V0(x,t)) 2 - m 2 2 (x,t)) - V(cf>(x,t)) 



+ E 2 



(«9o Xl (x, t)) 2 - (e- m Vxi(x, t)) 2 - /i 2 X 2 (x, t) - (x, t) 



12' 



Ey^(x,t)0 2 (x,t)}. 



(5.2) 



There exists an alternative derivation of the ensemble average of Eq. ( |3.7| ), which was 



presented for a single scalar field, as an intuitive argument in [|11|] . In the context of a single 
scalar field, the method is to work directly with the operator equation of motion for 0(x, t). 
The operator <p is reexpressed as the sum of a c-number tp c {t), representing the classical 
displacement, plus a shifted operator r/(x, t), 



(x,t) = <p e (t)+r}(x,t) 



(5.3) 



with ((j)(x,t))p = (p c (t). A thermal average is taken of this equation of motion, in which 
thermal expectation values involving ry(x, t) are computed such that (f c {t) is treated as an 
adiabatic parameter. To the order of perturbation theory considered in the previous sections, 
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for the single scalar field the intuitive derivation in ]TT| gives the same effective equation of 
motion as the ensemble average of Eq. ( |3. 7| ) as shown in || . 

No new considerations are needed to apply this intuitive derivation to the model in Eq. 
( p.ip . The x«( x , t) fields are treated as quantum fluctuations similar to ^(x, t). From the 
treatment in |H| , it follows that the expressions for m^, Ay and 7/1 in Eq. ( |3.7j ) will arise from 



the thermal averages, (?? 2 (x, t))p and (x 2 (x, t))p, taken with respect to the instantaneous 
background <p c {t). 

Although the approach in [|ll|] immediately isolates the dissipative term and the finite 



temperature renormalizations at the level of the equation of motion, it is not systematic to all 
orders. Furthermore, it cannot treat noise and it is valid only in the adiabatic approximation. 
These limitations can be accounted for in the closed-time-path formalism used in this paper. 
A recent work |36[ has discussed some of the difficulties associated with extending this 



formalism to an expanding background in order to treat noise and dissipation. Our goal at 
present is more modest. As an easier first step, the intuitive derivation of [JTTJ] is extended 
to an expanding background. 

The exact operator equations of motion from the Lagrangian Eq. ( |5.2j ) are 
0(x, t) + 3#0(x, t) - e" 2m V 2 0(x, t) + m 2 0(x, t) + — — — + £ <? 2 0(x, t) X 2 (x, t) = , 

(5.4) 

and 

Xi(x, t) + SHU*, t) - e- 2m V 2 Xi(x, t) + g 2 iXi {x, t)0 2 (x, t) + /i 2 Xi (x, t) + ^(x, t) 3 = . 

o 

(5.5) 

The objective is to displace the operator 0(x, t) by a x-independent c-number at time 
t = 0, (0(x, t = 0))p = (p c (0), and then determine the evolution of the expectation value 
(<j)(x.,t))p = (f c (t) by solving Eqs. ( |5.4| ) and (|5.5| ) perturbatively. Thus 0(x, t) is reexpressed 
as Eq. (|5.3j). With this definition of <f c {t), for flat, k — 0, nonexpanding, H = 0, spacetime, 
the resulting equation of motion is the same as the ensemble average of the equation of 



motion, Eq. (|3.71) . 
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For the case of expanding spacetime, H ^ 0, in order to obtain the equation of motion 
for <p c {t), thermal expectation values must be taken of Eqs. ( |5.4j ) and ( p.5| ). Provided the 
temperature, 1//3, of the thermal bath is time independent, (i.e., rapid equilibration time 
scales), thermal expectation values of terms linear in 0(x, t) can be replaced by (f c (t), just 
as for the nonexpanding case. In evaluating (r? 2 (x, t))p and (Xi( x ; if the characteristic 
time scale for the quantum fluctuations is much faster than the expansion time scale, l/H, 
the calculation is no different from the Minkowski space situation. This criteria is self- 
consistently satisfied provided 

r x ,rv»#, (5.6) 

where the left hand side is given in Eqs. ( [2.14|) and Q2.15] ) for our model. 



These arguments suggest that at leading nontrivial order, the effective equation of motion 
for ip c {t) in an expanding de Sitter spacetime, under the same conditions required for Eq. 
( p.7| ) plus the additional condition Eq. ( p.6| ) is 

<Pc(t) + {vivKt) + 3H)cp c {t) + m 2 TVc {t) + -£<p 3 e (t) = 0. (5.7) 

o 

Further justification that Eq. (|5.7|) is the appropriate replacement of Eq. (p.7|), for expanding 
de Sitter space, can be obtained from ||, where an effective equation of motion similar to 
Eq. (|5.7|) was obtained for a model like Eq. (|5.2f ). However, the coupling between fields was 
linear, (pXh which is analytically much more tractable than the present case of quadratic 
coupling, 4> 2 x 2 i- 

The entire discussion above assumes that the temperature has a well defined meaning in 
an expanding background. Furthermore, Eq. ( |5.7|) has been motivated under the restriction 
Eq. ( |5.6|) . As will be discussed next, condition (|5.6|) is a specific example of a general 
microscopic property argued in to be a necessary condition for warm inflation. As such, 
when (|5.7| ) is applied to the warm inflation scenario, condition ( |5.6| ) imposes no additional 
restriction. 

The warm inflation picture requires that an order parameter, in a strongly dissipative 
regime, slowly rolls down a potential, liberating vacuum energy into radiation energy, p r . The 
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nonisentropic expansion which underlies warm inflation imposes that the rate of radiation 
production is sufficient to compensate for red-shift losses due to cosmological expansion, 

H » H (5.8) 

Pt 

To give meaning to temperature, the newly liberated radiation must thermalize at a scale 
Trad which is faster than the expansion scale, 

r rad > H. (5.9) 

Minimally this requires an energy transfer rate from vacuum to radiation that is faster than 
the expansion rate, which in our model implies the condition (|5.6| ). Thus, Eq. ( |5.9|) is 
necessary to justify a temperature parameter T, which, combined with condition ( |5.6|) , are 
sufficient to justify the arguments leading to Eq. ( |5.7| ). To completely justify a temperature 
parameter for an expanding background spacetime, it is required studying the thermalization 
of the radiation, once it is liberated. General arguments, as well as specific calculations 
at high temperature, indicate that this rate is set by the temperature r rad ~ aT, 



for some appropriate, model dependent, coefficient a. This minimally requires T ^> H. 
However, a may be very small, as for example in Eqs. ( 2.14j) and (|2.15|) . Thus the correct 
constraint is 

aT > H. (5.10) 

This problem will not be considered further here. Eq. ( |5.6|) will be our only criteria for 
thermalization. This is equivalent to assuming that the thermalization rate is at least as 
fast as the energy transfer rate. 

Once Eq. ( |5.7| ) is accepted as the macroscopic equation governing the evolution of the 
order parameter <p c (t), it can be used as a given input to construct warm inflation scenarios 
as in J?|||. The microscopic origin of the equation can be forgotten up to restrictions on 
parameters and the self consistency condition Eq. ( [5.61 ). For a general equation like Eq. 
the warm inflation scenario requires the strong dissipative regime @: 
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[ri(<p e ) + 3H] ip c > (p c , (5.11) 

with i](ip c ) = rji(fc for our model. For the derivation in Sees. II-IV, where H = 0, this 
condition is sufficient to satisfy the adiabatic condition, Eq. (|4.2|) , which is required for the 
consistency of the microscopic calculation. As such, this model provides an example of a 
general point conveyed in ||, that warm inflation defines a good regime for application of 
finite temperature dissipative quantum field theory methods. The study of warm inflation in 
H|7].|J also found that to satisfy observational constraints on the expansion factor, it requires 

Tj((p c ) > 3H. (5.12) 

Thus, warm inflation is an extreme example of dissipative dynamic during de Sitter expan- 
sion. As demonstrated in PJ39|, dissipation is generally prevalent during inflation. The 



microscopic model in this paper could be used to examine the general case, but then the 
condition ( |5.12 ) can be relaxed. Here, only the warm inflation regime will be further exam- 



ined. Thus in the limit given by Eq. ( |5.12| ) and based on the remaining discussion in this 
section, the equation of motion for the order parameter ip c (t) in our model for the warm 
inflation scenario turns out to be Eq. (|4.4j) but with the additional constraint Eq. (|5.6|) . 

The other input for constructing warm inflation scenarios is the free energy in the ex- 
panding environment for the model (|5.2|) . It already has been argued above that temperature 
is a good parameter for describing the state of the radiation in the warm inflation regime. 
It also follows from the above that the change in temperature can be treated adiabatically 
in the thermodynamic functions, since this requires r ra( j ^> T/T, which is automatically 
satisfied due to Eq. ( |5.9| ). Therefore the free energy density should be well represented by 
the Minkowski space expression, with temperature treated as an adiabatic parameter. For 
the model in Sees. II-IV, the free energy density is 

jpf rri\ m T 2 i 4 (N+ 1)tt 2 4 

F{<Pc T) = — Vc + — Vc — T , (5.13) 

where the factor N + 1, in the last term, comes from the functional integration over the \ 
fields and the 0-field's fluctuations. Having established this to be the free energy density 
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for the warm inflation scenario, the other thermodynamic functions such as pressure, energy 
density and entropy density can be easily obtained. 

With the free energy ( |5.13| ) and the order parameter equation of motion, Eq. ( p.7| ), 



determined, the time evolution of the three unknowns: temperature T(t), scale factor R(t) 
and order parameter ip c (t), can be obtained from Eq. ( p-7| ) plus any two independent 
equations from Friedmann cosmology along with a self consistency check for adiabaticity, Eq. 
( |5.11| ). At this point the procedure in || can be followed. However due to the microscopic 



origin of this model, additional self consistency checks are necessary for adiabaticity, given 
by Eq. ( |4.2|) and thermalization, Eq. (|5.6| ). Observationally interesting expansion factors 
will require H > 0/0, in which case the condition (|5.6|) immediately implies the microscopic 



adiabatic condition (|4.2|). 



B. Results 

Up to this point, the formulation of warm inflation in conjunction with a microscopic 
dynamics has been general. In the remainder of this section, some demonstrative calculations 
of this cosmology will be presented based on our microscopic model. An exhaustive analysis 
of the parameter space will not be performed. In this first examination, the emphasis is to 
understand the interplay between the microscopic and macroscopic physics of warm inflation 
for generic potentials, which in particular, have curvature scale of order the temperature 
scale. For such potentials, thermal fluctuations that displace y2 c (0) substantially from the 
origin are exponentially suppressed. However, it is such fluctuations that allow enough 
time, during the roll down back to the origin, for the universe to inflate sufficiently As such, 
this elementary fact, in any case, quells significant interest in comparing the cases we will 
examine to observation. 

It should be noted that the order parameter in this symmetry restored warm inflation 



regime is configured similar to those in the chaotic inflation scenario ||40|| . However, in the 



chaotic inflation scenario the potentials are ultra-flat. Such potentials permit large fluctu- 
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ations of the order parameter and in fact prefer them. The dissipative model in this paper 
could be studied for the case of ultra-flat potentials, perhaps motivated by supersymmetric 
model building. This would extend the pure quantum mechanical, new-inflation type dy- 



namics of chaotic inflation into the intermediate regime discussed in P,p9||. This will not be 
examined here. 

Proceeding with our demonstrative examination of warm inflation, let the origin of time 
be the beginning of the inflationlike regime (BI) and also the beginning of our treatment. 
The basic picture of the particular warm inflation scenario studied here is as follows. At 
t = the initial conditions are arranged so that the field is displaced from the origin 
(0(0)) = 4>Bii the temperature of the universe is T BI and since the universe is at the onset 
of the inflation-like regime, by definition this means the vacuum energy density equals the 
radiation energy density, p v (0) = p r (0). For t > the field will relax back to the origin 
within a strongly dissipative regime and in the process liberate vacuum energy into radiation 
energy. Simultaneously, the scale factor will undergo inflation- like expansion. During the 
roll-down period, the vacuum energy first dominates until at some point it is superseded by 
the radiation energy. At this point the universe smoothly exits the inflation- like regime into 
the radiation dominated regime. 

From our model in the previous sections, we will consider the case of N' x-bosons (x') 
with gj = g > fa, j = 1 . . . N' and N - N' x-bosons (x) with gj <^fj = f,j = N' + l...N. 
For this model, the dissipative dynamics of <p c , expressed through m T , X T and T]i, is controlled 
by the N' x-bosons. The other N—N' x-bosons only serve as additional fields in the radiation 



bath. For this purpose, from Eq. ( 2.15|) , for / > g 2 the x an d x' bosons will be equally 



effective in thermalizing the radiation energy. 

In this paper we will examine this scenario in the regime 

» (5.M) 

and in the high temperature limit T ^> my, /iy. Also, for ease of presentation, we will write 
the expression for T(q) at q = 0. Although with these simplifications the results will not be 
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cosmologically interesting, it is a good example to demonstrate the general procedure. In 
this regime, the effective equation of motion for (p c , from Eq. ( [4.4| ), is: 



dt 4 

where 



-ip c , (5.15) 



2 = 3^7 ~ 72NWZr B i/f>r Bi y (5 ' 16) 
Formally the Friedmann cosmology for the warm inflation scenario associated with the above 
equation was called the quadratic limit in |§. 

The macroscopic and microscopic requirements of warm inflation will imply various para- 
metric constraints which are as follows. Eq. ( |5.14j ) will be satisfied by requiring: 

24 2 y ' ' 

where the parameter r> 1 has been introduced. As shown in ||, throughout the inflation- 
like period until just before it ends, the temperature drops slightly faster than (f>. As such, 
the thermal mass term, m^0 2 /2 ~ T 2 <ft 2 , will continue to satisfy Eq. (|5.14|) given that 
initially it does. 

The number of e- folds, N e obtained during the roll-down is, from ||, 

Ne*^, (5.18) 

where 



H 



8tt\ t <P bi 



(5.19) 



fr_ >> {5 20) 

1927T 3m p 



\ 72m2 ' 

with m p the Planck mass. The microscopic condition, Eq. ( |5.6|) , requires 



The threshold condition for inflation, p v (0) = p r (0), implies 



' JV+ ^V B; = ^k. (5.21) 
30 BI 24 v ; 
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Finally, the validity of the perturbative derivation in the previous sections will require 

g,X^,X T <\, (5.22) 

where m 2 T ~ (A + 2N'g 2 )T 2 /2A and X T ~ A - ^ff-N'g 3 . 

Eq. Q5.20D can be turned into an equality, in which case, along with Eqs. ( pM7| ), fl5.18| ). 



and C |5 . 2 1|) , they determine the boundary of the allowed parameter space. Thus there are six 
constraining equations for the eleven quantities A, g, At, mr, m p , <Pbi, N, N', N e , r, and Tbi- 
We will let Tbi set the overall scale and will fix N' , N e , r, g. Then, based on the constraint 
equations, this determines the remaining parameters. In particular, we have for At: 

A , 3N'g*\n(2y/U/g) 



4iV> 2 



(5.23) 



This expression is suggestive of the case 3 analyzed in Sec. IV. For such, taking then 



(3gy/3/2ir + 2^ N' g 2 T 2 /24, we get the additional parameters: 

2nT BI 



{3gVS/2n + 2) N e i 
g \| 6 ln(2VT2/0) 



(5.24) 



Ar A 5r 2 N'N e f3V3 \ 

N+l^ — \-—g + 2\ , 5.25 

12 \n(2y/U/g) \ 2vr y J K J 



and 



192vr 2 rT B/ (3y/3 , 



3nN'N e 



(5.26) 



9g 4 \ 2tt y J \\ \n{2Vl2/g) ' 

Based on these equations, it is not difficult to find parametric regimes in which the warm 
inflation scenario is realized, but it is only for N e < 1. As such, this simple case has no 
observational relevance. There are a few improvements that could be made to our analysis 
that would increase N e . Firstly our estimates above ignore the effects of the thermal mass 
term, mf^/2, on the dynamics. Its contribution to the energy density and pressure are 
Pm T = — ^t^c/2 an d Pm T = Pm T respectively. Thus it helps the term to drive inflation. 
Secondly, recall that the dissipative coefficient rjx ~ 1/T. In the above analysis, we fixed 
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T = Tbi- However, during the roll-down, temperature does fall by a factor of order ten, 
which in turn would increase 771. Finally the parameter regime could be extended to include 
both positive and negative 4> — x couplings, such as the ASR case described in subsection 
IV.A. As noted there, in this regime the duration of overdamped motion can be increased 
significantly within the perturbative amplitude expansion. This directly corresponds to 
increasing the e- folds N e . 

A more elementary modification is to extend the region of validity to larger displacements 
of ip c . The extension to this larger regime can be treated by a summation of the complete one- 
loop series as outlined in subsection IV. B. Overdamped motion for much larger displacements 
of ip c can also be attained by a modification to the <p — x interaction in Eq. (|2.1| ) as 



In this distributed mass model (DMM), a given Xj-field is thermally excited when its effective 
mass g 2 (ip c — Mj) 2 < T 2 . The contribution from the thermally excited x-Aelds to the effective 
dynamics of tp c is similar to our calculations in sections II-IV. As such, an effective equation 
of motion similar to Eq. (4.4) can be obtained for this modified model. Given an appropriate 
distribution of mass coefficients Mj along the path of ip c , from an arbitrarily large initial 
displacement ifBi, <f c {t) could undergo overdamped motion along its entire path. Details 
will be presented elsewhere on the warm inflation scenario which considers these various 
cases. 

For this "symmetry restored" case, initial fluctuations of (p c (t), ¥>bi, are strongly sup- 
pressed with probability exp(— volume p/Tbi) ~ exp[— (128t[) 3 N N e r 2 / g 12 }, where the volume 
~ 1/H 3 . The most optimistic initial conditions have probability ~ exp[— 1 x 10 7 ]. Thus, 
unless a viable mechanism is found to justify a large enough initial value of <p c , the regime 
investigated here may not be very relevant for practical applications of warm inflation. In 
any case, the microscopic dynamics of the symmetry restored regime investigated here is 
similar to more realistic scenarios in the symmetry broken regime, where the field has an 
average initial value close or identical to zero. An important difference being that the initial 




(5.27) 
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state in the latter case has no Boltzmann suppression. 

This section has made an initial examination of treating strong dissipation from first 
principles during a de Sitter expansion regime. Further results will be presented elsewhere 
as well as a calculation similar to this one, for the symmetry broken case. 



VI. CONCLUSION 

In this paper a microscopic quantum field theory model has been presented, describing 
overdamped motion of a scalar field. Commonly, such behavior is treated phenomenolog- 
ically by Ginzburg-Landau order parameter kinetics. Our model provides a first principle 
explanation of how kinetics equivalent to the Ginzburg-Landau type, which is first order in 
time, arise for inherently second order dynamical systems. The microscopic treatment of 
this problem, in principle, should be well controlled, due to its fundamental reliance on the 
adiabatic limit, and our model exemplifies this expectation. 

The calculational method for treating dissipation in this paper has one distinct difference 
from several other related works. In our calculation, we consider the effect of particle 
lifetimes in the effective equation of motion. To our knowledge, this effect has been discussed 
in only a few works in the past [p|, ll|JH • 



Sees. II-IV presented a general, flat-space treatment, which offers a microscopic jus- 
tification to the often used limit of diffusive Ginzburg-Landau scalar field dynamics. We 
have shown how it is possible to obtain an effective evolution for the scalar field which is 
first-order in time, due to its own thermal dissipation effects, interpreted microscopically 
as its decay into many quanta. In a sense, the field acts as its own brakes, the slowing of 
its dynamics being attributed to the highly viscous medium where it propagates, a densely 
populated sea of its own decay products. 

The application that we considered in Sec. V was in expanding spacetime, for the 
cosmological warm inflation scenario. Although we did not formally derive the extension of 
our flat-space model of Sees. II-IV to an expanding background, we did present heuristic 
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arguments that validate this extension for the special needs of warm inflation. The results 
of the simple analysis in Sec. V are strongly dependent on initial conditions and may be 
difficult to implement for models of observational interest. Nevertheless, these results will 
provide useful guidance both for modifications of this model and for our next study of the 
symmetry broken regime. 

The direct significance of the present study to inflationary cosmology would be to the 



initial state problem jy] in scenarios during symmetry breaking. The initial conditions 
required for warm inflation in the symmetry broken case are similar to new inflation. The 
requirement is a thermalized inflaton field, which at the onset of the warm inflation regime 
is homogeneous with expectation value ((f)) p = 0. Although we have made no detailed 
application of our results to this problem, some general features are evident from the analysis 
in Sec. V. In particular, both the suppression of large fluctuations and thermalization are 
mutually consistent with strong dissipative dynamics. Many of the difficulties that have 
been discussed [^1,^2] in association with the initial state problem, are eliminated in the 
strong dissipative regime. In addition, the damping of fluctuations should simplify the 
formal problem of coupling this model to classical gravity. Thus, the strong dissipative 
regime appears to have the correct features both to carry the universe into an inflation-like 
phase and then to smoothly exit into a hot big-bang regime. 
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APPENDIX A: 



We now give a brief overview of the calculation of the imaginary part of the two-loop, 
setting sun, self-energy terms in ( |2.10| ) and ( |2.11| ). Let us express generically those diagrams 
in terms of field propagators with masses m s and m t and external lines of type s. For 
an interaction between s and t fields of the form g 2 t /2 s 2 t 2 (for s — t — <p {Xj), 9st = 



A/12 (A/12) and for s 



Xj, g 2 st = gf) the imaginary part of the two- loop sunset 



diagram for the s field can be expressed by (see for example, in Ref. ||21|| , Appendix G and 



also Wang and Heinz in |22|) 
ImS s (g) = Im [ _ 




d 3 kd 3 l 



a=±l 



rrt e s E l 

^k+q^k+l 



1 + n(a k+q E s k+q ) 



X 



1 + nfaEt) 1 + n(a k+l E{ +l ) 6(E s q - a x E\ - a k+q E s k+q - a k+l E t 



4 



(Al) 



where Ef, = yk 2 + m 2 , n(E) is the bose-distribution function and S s j is a symmetry factor: 
for s = t, S = 12 and for s ^ t, S = 1. By expanding in the sum in a and retaining only 
the on-shell, energy conserving processes (corresponding to the scattering processes st — > st, 
ts — ► ts), we obtain the result: 

d 3 kd 3 l 



ImS s (g) = S„ t tgi t v 



1-e 



4(2tt) 5 

x 5(E s q - E\ - E s k+q + El +l ) 



1 + n(E s k+q ) 1 + n{E\) n(E 



A+i) 



(A2) 



with S Sj t = 1 for m s = m t and 5 s j = for m s ^ m t . ImS(g) for s = t has been obtained 



in details in Refs. pTJ and ||22|| . In particular, Wand and Heinz in [^2j have discussed 
and obtained in detail the kinematic limits, for the A0 4 model, of the integration on the 
momenta in (|A2|), implicit in the Dirac delta-function. Here we obtain the results for the 
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case of m s ^ m t . In (|A2|) , by defining in the three-dimensional momentum integrations the 
angular differentials d(cos9k) = dE^ + q E% + J (kq) and d(cos9i) = dEl +l El +l /(kl) (k,q,l = 
|k|, |q|, |1|), we are then able to perform the angular integrals in ( |A2| ) . From the constraint 
in the integration limits for k and I, which comes from the delta-function, we obtain the 
following result for ImE s (g), 



Im£ s (g) = S S;t gt t 



T (l + ±5 S 



(1 



x In 



4(2?r) 3 g 
ePH+i - 1 



dk 



o 



eP&k+i - e m- E i) 

1 + n{E s q - Ef)} \l + n(El 



«(-*.«) Idl 

—u(k,q) Ei 
Idl 



1 + n( E s q - Eft] \l + n{E\ 



IdV 



OO roc 

dk / ^ — I dk . 

-q J u(k+q,q) -C/£ Jq Ju(—k+q,q) -C/; 



In 



- 1 



J3{E'-Ef) 



(A3) 



where, in the above expressions, the function u(k, q) is given by 

1 



u(k,q) 



2 ^ kn> 



(k — q) 2 + m 2 s — Jq 2 + m 2 



(k 2 m 2 s + 2m 2 s m 2 + 2q 2 m 2 



2kqm t + 2m t J (k — q) 2 + m 2 Jq 2 + m 



(A4) 



In the first term of ( |A3| ) , we can make the change of integration variables k + l — k', I — I', to 
obtain VdV /E\, = dEj,. Doing the same for the remaining integrals, Idl/Ej = dEf, relabeling 
I' and k' back to /, k, and performing a new change of integration variables, y = e~@ E i , 
dy = -Pe- 0E idEf, we are able to compute the y integrations. From this point on, the 



integrations are equivalent to the ones in [21], once we change the integration limits in the 



k and / integrals, and take into account the function u(k,q), Eq. (|A4|) . The results shown 
in Sec. II for r^, r x ., Eqs. (|2.14j ) and Q2.15| ), are obtained once the limit of vanishing q 
(|q| = 0), is taken in the above equations. 



APPENDIX B: 



Let us now compute the expression on the LHS of ( |3.5| ) in the adiabatic approximation, 
stressing the need for regularizing the propagators with quasi-particle lifetimes. Let us write 
the LHS of ( |3.5|) in terms of the adiabatic series: 
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[ d A x'ip 2 c (x')Im \G ++ } 2 0(t-t') 

J - - 1 x,x' 

= / d~x 



n 



(2vr) 3 n! <9t 



/ 71 



, =i y (2tt) 



d3g ^ dt\t'-t) n 



oo 



Im[G ++ (q, |t-t'|)G ++ (q-k, |t - , (Bl) 



where G ++ (q, t — t') can be read from ( |2.20| ) and ( |2.21| ) and it refers generically to the <fi or 



Xj field propagators. The approximation of considering a homogeneous field, (p c = <p c {t), is 
equivalent to taking the limit k — > for the external momentum in G ++ (q— k, \t—t'\). In this 
case the x' and k integrations in ( [BID can be done trivially. However, it is known that taking 
the limit of zero external momentum (/c M — > 0) in self-energy expressions, which are 

given in terms of products of non-local propagators, can be problematic. This is related to 
the non-analyticity of these expressions around the origin. In particular, Gross, Pisarski, 
and Yaffe in |^3[ argue that a correct way of taking the limit k^ — > 0, at finite temperatures, 



is to first take ko — > and then k — > 0. They also argue that the non-analyticity problem, 
which comes from a failure to do a self- consistent calculation, would be eliminated once 
fully-dressed propagators are taken consistently into account. We note, in particular, that 
for fully-dressed propagators the decay width T works Kulator. 

In ( |B1[ ), the limit k — > is implicit in the adiabatic approximation, where the fields are 
required to change slowly in time, while the limit k — > is implemented by approximating 
the fields to be homogeneous. In evaluating ( |B1|) , we will first compute the time integral, 
expand in terms of the "regulator" T, and then finally take the limit k — > 0. This is 
the opposite of what was done in obtaining Eq. (|3.5| ). From this lesson we will see both 
the importance of considering fully-dressed propagators and how consistent results can be 
obtained once the limits are taken in the correct order. This will be crucial when evaluating 
the dissipation contribution (n — 1, in ([Bl])). 

Consider the time integral on the RHS of Eq. ( |B1| ) and expand it to first order in the 
"regulator", T (which is of order A 2 , g 4 ). Higher order terms in T need to be considered in 
conjunction with higher order loop terms, for consistency. For n = 0, the zeroth-order in 
the adiabatic approximation, we obtain 
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dt'Im [G ++ (q, \t - t'\)G ++ (q -k,\t- t'\) 



+ ^q-k 



+ 2 



w q -ki(w q ) - w q n(w q _ k ) 



4w q u; q „ k 

Taking the limit k — ► 0, we obtain 



^q - W q-k 



+ 



r Q + r 



q-k 



s ^q + ^q-k, 



(B2) 



lim / dt'Im 

k^0.7-OC 



G++(q,|t-t'|)G++(q-k,|t-i'|) 



[1 + 2n(w q )] _ ^ n(a; q ) [1 + w(q; q )] + Q 



50^ 



q 



(B3) 



As expected, ( p3|) is recognized to be the usual 1-loop correction to the quartic coupling 
constant. 

For n=l, the first order in the adiabatic approximation, which will give the dissipation 
coefficient, we obtain (again retaining terms up to first order in T) 

dt'{t' - t)Im [G ++ (q, \t - f|)G ++ (q - k, \t - t'\) 



[1 + 2n(cu q )] 



3^ q + ^q-k 



Tq + Tq-k 



W q K - W q-k 



) K-^ q _k)^ + ^- 



+ (u; q ^ w q _ k J 



-/5r q n(^ q ) [1 + ra(o> q )] 



(Ua # ^q-k) + 



r 



q-k 



W q + Wq-k. 



(B4) 



We can see, in contrast with (|B2|) , that the limit k — > is divergent. This divergence is 
reminiscent of the on-shell singularity which is present in the integral in ([B4f) when bare 
propagators are used, thus showing the importance of T as a regulator. By first taking the 
homogeneous limit k — > and then expanding in T, we obtain the result given in the text, 
which is the first term on the RHS in ( |3.5|) . 

The same calculation can be performed for the n = 2 case, the second order in the 
adiabatic expansion, from which we obtain 

dt'{t' - t) 2 Im [G ++ (q, \t - t'|)G ++ (q - k, \t - t'\) 



lim 

k^O j-oo 



i + 2nK) | c /T 2 



16u; q 



(B5) 



This result is consistent with a recent calculation in Ref. ML which addresses the time 



derivative expansion of the effective action, for a given scalar field model. The authors of 
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|44| also discuss their work in the context of the non-analyticity problem in finite temperature 



QFT. In fact, we note that the second order term in the adiabatic approximation can be 
associated with the first order term in the time derivative expansion of the effective action, 
~ Z(ip)(d t (p) 2 (for a time- dependent, space-homogeneous field configuration). 
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FIGURE CAPTIONS 



Figure 1: ImE x (q, a; x (q)) normalized by its |q| = value, for different values of masses 
and space momentum. 

Figure 2: The dissipation coefficient rji computed (for N = 25) with ImE x (q, u; x (q)) 
and with the approximation |q| = for the imaginary part of the self-energy. 

Figure 3: Results for the adiabatic condition, Eq. ( [4.71 ). The dashed lines correspond 
to \(p c /(p c \, for different values for (p c . The solid line corresponds to T x (g), evaluated at 
|q| = T. All cases shown are for N = 25. The region satisfying the adiabatic condition is 
the intersection of the region above the dashed lines with the region below the solid line. 
Fig. 3a is for A = g 2 , Fig. 3b is for A = g and Fig. 3c is for X T = g A . 
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Figure 2 
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Figure 3b 
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Figure 3c 
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